An Adult Zebrafish Model Reveals that Mucormycosis Induces Apoptosis of Infected Macrophages.

Mucormycosis is a life-threatening fungal infection caused by various ubiquitous filamentous fungi of the Mucorales order, although Rhizopus spp. and Mucor spp. are the most prevalent causal agents. The limited therapeutic options available together with a rapid progression of the infection and a difficult early diagnosis produce high mortality. Here, we developed an adult zebrafish model of Mucor circinelloides infection which allowed us to confirm the link between sporangiospore size and virulence. Transcriptomic studies revealed a local, strong inflammatory response of the host elicited after sporangiospore germination and mycelial tissue invasion, while avirulent and UV-killed sporangiospores failed to induce inflammation and were rapidly cleared. Of the 857 genes modulated by the infection, those encoding cytokines, complement factors, peptidoglycan recognition proteins, and iron acquisition are particularly interesting. Furthermore, neutrophils and macrophages were similarly recruited independently of sporangiospore virulence and viability, which results in a robust depletion of both cell types in the hematopoietic compartment. Strikingly, our model also reveals for the first time the ability of mucormycosis to induce the apoptosis of recruited macrophages but not neutrophils. The induction of macrophage apoptosis, therefore, might represent a key virulence mechanism of these fungal pathogens, providing novel targets for therapeutic intervention in this lethal infection.

Molecular Tools for Carotenogenesis Analysis in the Mucoral Mucor circinelloides.

The carotene producer Mucor circinelloides is the fungus within the Mucoromycota phylum with the widest repertoire of molecular tools to manipulate its genome. The initial development of an effective procedure for genetic transformation and later improvements have resulted in an expansion of available tools, which include gene replacement, inactivation of gene expression by RNA silencing, gene overexpression, and functional genomics. Moreover, sequencing of its genome has given a definitive boost to these techniques making attainable the study of genes involved in many physiological or developmental processes, including carotenoid biosynthesis. Here, we describe in detail the latest molecular techniques currently used in M. circinelloides that have made it a valuable model for studying gene function within its phylum.

The RNAi Universe in Fungi: A Varied Landscape of Small RNAs and Biological Functions.

RNA interference (RNAi) is a conserved eukaryotic mechanism that uses small RNA molecules to suppress gene expression through sequence-specific messenger RNA degradation, translational repression, or transcriptional inhibition. In filamentous fungi, the protective function of RNAi in the maintenance of genome integrity is well known. However, knowledge of the regulatory role of RNAi in fungi has had to wait until the recent identification of different endogenous small RNA classes, which are generated by distinct RNAi pathways. In addition, RNAi research on new fungal models has uncovered the role of small RNAs and RNAi pathways in the regulation of diverse biological functions. In this review, we give an up-to-date overview of the different classes of small RNAs and RNAi pathways in fungi and their roles in the defense of genome integrity and regulation of fungal physiology and development, as well as in the interaction of fungi with biotic and abiotic environments.

A non-canonical RNA degradation pathway suppresses RNAi-dependent epimutations in the human fungal pathogen Mucor circinelloides.

Mucorales are a group of basal fungi that includes the casual agents of the human emerging disease mucormycosis. Recent studies revealed that these pathogens activate an RNAi-based pathway to rapidly generate drug-resistant epimutant strains when exposed to stressful compounds such as the antifungal drug FK506. To elucidate the molecular mechanism of this epimutation pathway, we performed a genetic analysis in Mucor circinelloides that revealed an inhibitory role for the non-canonical RdRP-dependent Dicer-independent silencing pathway, which is an RNAi-based mechanism involved in mRNA degradation that was recently identified. Thus, mutations that specifically block the mRNA degradation pathway, such as those in the genes r3b2 and rdrp3, enhance the production of drug resistant epimutants, similar to the phenotype previously described for mutation of the gene rdrp1. Our genetic analysis also revealed two new specific components of the epimutation pathway related to the quelling induced protein (qip) and a Sad-3-like helicase (rnhA), as mutations in these genes prevented formation of drug-resistant epimutants. Remarkably, drug-resistant epimutant production was notably increased in M. circinelloides f. circinelloides isolates from humans or other animal hosts. The host-pathogen interaction could be a stressful environment in which the phenotypic plasticity provided by the epimutant pathway might provide an advantage for these strains. These results evoke a model whereby balanced regulation of two different RNAi pathways is determined by the activation of the RNAi-dependent epimutant pathway under stress conditions, or its repression when the regular maintenance of the mRNA degradation pathway operates under non-stress conditions.

RNAi-Based Functional Genomics Identifies New Virulence Determinants in Mucormycosis.

Mucorales are an emerging group of human pathogens that are responsible for the lethal disease mucormycosis. Unfortunately, functional studies on the genetic factors behind the virulence of these organisms are hampered by their limited genetic tractability, since they are reluctant to classical genetic tools like transposable elements or gene mapping. Here, we describe an RNAi-based functional genomic platform that allows the identification of new virulence factors through a forward genetic approach firstly described in Mucorales. This platform contains a whole-genome collection of Mucor circinelloides silenced transformants that presented a broad assortment of phenotypes related to the main physiological processes in fungi, including virulence, hyphae morphology, mycelial and yeast growth, carotenogenesis and asexual sporulation. Selection of transformants with reduced virulence allowed the identification of mcplD, which encodes a Phospholipase D, and mcmyo5, encoding a probably essential cargo transporter of the Myosin V family, as required for a fully virulent phenotype of M. circinelloides. Knock-out mutants for those genes showed reduced virulence in both Galleria mellonella and Mus musculus models, probably due to a delayed germination and polarized growth within macrophages. This study provides a robust approach to study virulence in Mucorales and as a proof of concept identified new virulence determinants in M. circinelloides that could represent promising targets for future antifungal therapies.

A new regulatory mechanism controlling carotenogenesis in the fungus Mucor circinelloides as a target to generate ß-carotene over-producing strains by genetic engineering.

BACKGROUND: Carotenoids are natural pigments with antioxidant properties that have important functions in human physiology and must be supplied through the diet. They also have important industrial applications as food colourants, animal feed additives and nutraceuticals. Some of them, such as ß-carotene, are produced on an industrial scale with the use of microorganisms, including fungi. The mucoral Blakeslea trispora is used by the industry to produce ß-carotene, although optimisation of production by molecular genetic engineering is unfeasible. However, the phylogenetically closely related Mucor circinelloides, which is also able to accumulate ß-carotene, possesses a vast collection of genetic tools with which to manipulate its genome. RESULTS: This work combines classical forward and modern reverse genetic techniques to deepen the regulation of carotenoid synthesis and generate candidate strains for biotechnological production of ß-carotene. Mutagenesis followed by screening for mutants with altered colour in the dark and/or in light led to the isolation of 26 mutants that, together with eight previously isolated mutants, have been analysed in this work. Although most of the mutants harboured mutations in known structural and regulatory carotenogenic genes, eight of them lacked mutations in those genes. Whole-genome sequencing of six of these strains revealed the presence of many mutations throughout their genomes, which makes identification of the mutation that produced the phenotype difficult. However, deletion of the crgA gene, a well-known repressor of carotenoid biosynthesis in M. circinelloides, in two mutants (MU206 and MU218) with high levels of ß-carotene resulted in a further increase in ß-carotene content to differing extents with respect to the crgA single-null strain; in particular, one strain derived from MU218 was able to accumulate up to 4 mg/g of ß-carotene. The additive effect of crgA deletion and the mutations present in MU218 suggests the existence of a previously unknown regulatory mechanism that represses carotenoid biosynthesis independently and in parallel to crgA. CONCLUSIONS: The use of a mucoral model such as M. circinelloides can allow the identification of the regulatory mechanisms that control carotenoid biosynthesis, which can then be manipulated to generate tailored strains of biotechnological interest. Mutants in the repressor crgA and in the newly identified regulatory mechanism generated in this work accumulate high levels of ß-carotene and are candidates for further improvements in biotechnological ß-carotene production.

RNAi pathways in Mucor: A tale of proteins, small RNAs and functional diversity.

The existence of an RNA-mediated silencing mechanism in the opportunistic fungal pathogen Mucor circinelloides was first described in the early 2000. Since then, Mucor has reached an outstanding position within the fungal kingdom as a model system to achieve a deeper understanding of regulation of endogenous functions by the RNA interference (RNAi) machinery. M. circinelloides combines diverse components of its RNAi machinery to carry out functions not only limited to the defense against invasive nucleic acids, but also to regulate expression of its own genes by producing different classes of endogenous small RNA molecules (esRNAs). The recent discovery of a novel RNase that participates in a new RNA degradation pathway adds more elements to the gene silencing-mediated regulation. This review focuses on esRNAs in M. circinelloides, the different pathways involved in their biogenesis, and their roles in regulating specific physiological and developmental processes in response to environmental signals, highlighting the complexity of silencing-mediated regulation in fungi.

Distinct RNAi Pathways in the Regulation of Physiology and Development in the Fungus Mucor circinelloides.

The basal fungus Mucor circinelloides has become, in recent years, a valuable model to study RNA-mediated gene silencing or RNA interference (RNAi). Serendipitously discovered in the late 1900s, the gene silencing in M. circinelloides is a landscape of consensus and dissents. Although similar to other classical fungal models in the basic design of the essential machinery that is responsible for silencing of gene expression, the existence of small RNA molecules of different sizes generated during this process and the presence of a mechanism that amplifies the silencing signal, give it a unique identity. In addition, M. circinelloides combines the components of RNAi machinery to carry out functions that not only limit themselves to the defense against foreign genetic material, but it uses some of these elements to regulate the expression of its own genes. Thus, different combinations of RNAi elements produce distinct classes of endogenous small RNAs (esRNAs) that regulate different physiological and developmental processes in response to environmental signals. The recent discovery of a new RNAi pathway involved in the specific degradation of endogenous mRNAs, using a novel RNase protein, adds one more element to the exciting puzzle of the gene silencing in M. circinelloides, in addition to providing hints about the evolutionary origin of the RNAi mechanism.

The RNAi machinery controls distinct responses to environmental signals in the basal fungus Mucor circinelloides.

BACKGROUND: RNA interference (RNAi) is a conserved mechanism of genome defence that can also have a role in the regulation of endogenous functions through endogenous small RNAs (esRNAs). In fungi, knowledge of the functions regulated by esRNAs has been hampered by lack of clear phenotypes in most mutants affected in the RNAi machinery. Mutants of Mucor circinelloides affected in RNAi genes show defects in physiological and developmental processes, thus making Mucor an outstanding fungal model for studying endogenous functions regulated by RNAi. Some classes of Mucor esRNAs map to exons (ex-siRNAs) and regulate expression of the genes from which they derive. To have a broad picture of genes regulated by the silencing machinery during vegetative growth, we have sequenced and compared the mRNA profiles of mutants in the main RNAi genes by using RNA-seq. In addition, we have achieved a more complete phenotypic characterization of silencing mutants. RESULTS: Deletion of any main RNAi gene provoked a deep impact in mRNA accumulation at exponential and stationary growth. Genes showing increased mRNA levels, as expected for direct ex-siRNAs targets, but also genes with decreased expression were detected, suggesting that, most probably, the initial ex-siRNA targets regulate the expression of other genes, which can be up- or down-regulated. Expression of 50% of the genes was dependent on more than one RNAi gene in agreement with the existence of several classes of ex-siRNAs produced by different combinations of RNAi proteins. These combinations of proteins have also been involved in the regulation of different cellular processes. Besides genes regulated by the canonical RNAi pathway, this analysis identified processes, such as growth at low pH and sexual interaction that are regulated by a dicer-independent non-canonical RNAi pathway. CONCLUSION: This work shows that the RNAi pathways play a relevant role in the regulation of a significant number of endogenous genes in M. circinelloides during exponential and stationary growth phases and opens up an important avenue for in-depth study of genes involved in the regulation of physiological and developmental processes in this fungal model.

A non-canonical RNA silencing pathway promotes mRNA degradation in basal Fungi.

The increasing knowledge on the functional relevance of endogenous small RNAs (esRNAs) as riboregulators has stimulated the identification and characterization of these molecules in numerous eukaryotes. In the basal fungus Mucor circinelloides, an emerging opportunistic human pathogen, esRNAs that regulate the expression of many protein coding genes have been described. These esRNAs share common machinery for their biogenesis consisting of an RNase III endonuclease Dicer, a single Argonaute protein and two RNA-dependent RNA polymerases. We show in this study that, besides participating in this canonical dicer-dependent RNA interference (RNAi) pathway, the rdrp genes are involved in a novel dicer-independent degradation process of endogenous mRNAs. The analysis of esRNAs accumulated in wild type and silencing mutants demonstrates that this new rdrp-dependent dicer-independent regulatory pathway, which does not produce sRNA molecules of discrete sizes, controls the expression of target genes promoting the specific degradation of mRNAs by a previously unknown RNase. This pathway mainly regulates conserved genes involved in metabolism and cellular processes and signaling, such as those required for heme biosynthesis, and controls responses to specific environmental signals. Searching the Mucor genome for candidate RNases to participate in this pathway, and functional analysis of the corresponding knockout mutants, identified a new protein, R3B2. This RNase III-like protein presents unique domain architecture, it is specifically found in basal fungi and, besides its relevant role in the rdrp-dependent dicer-independent pathway, it is also involved in the canonical dicer-dependent RNAi pathway, highlighting its crucial role in the biogenesis and function of regulatory esRNAs. The involvement of RdRPs in RNA degradation could represent the first evolutionary step towards the development of an RNAi mechanism and constitutes a genetic link between mRNA degradation and post-transcriptional gene silencing.

Antifungal drug resistance evoked via RNAi-dependent epimutations.

Microorganisms evolve via a range of mechanisms that may include or involve sexual/parasexual reproduction, mutators, aneuploidy, Hsp90 and even prions. Mechanisms that may seem detrimental can be repurposed to generate diversity. Here we show that the human fungal pathogen Mucor circinelloides develops spontaneous resistance to the antifungal drug FK506 (tacrolimus) via two distinct mechanisms. One involves Mendelian mutations that confer stable drug resistance; the other occurs via an epigenetic RNA interference (RNAi)-mediated pathway resulting in unstable drug resistance. The peptidylprolyl isomerase FKBP12 interacts with FK506 forming a complex that inhibits the protein phosphatase calcineurin. Calcineurin inhibition by FK506 blocks M. circinelloides transition to hyphae and enforces yeast growth. Mutations in the fkbA gene encoding FKBP12 or the calcineurin cnbR or cnaA genes confer FK506 resistance and restore hyphal growth. In parallel, RNAi is spontaneously triggered to silence the fkbA gene, giving rise to drug-resistant epimutants. FK506-resistant epimutants readily reverted to the drug-sensitive wild-type phenotype when grown without exposure to the drug. The establishment of these epimutants is accompanied by generation of abundant fkbA small RNAs and requires the RNAi pathway as well as other factors that constrain or reverse the epimutant state. Silencing involves the generation of a double-stranded RNA trigger intermediate using the fkbA mature mRNA as a template to produce antisense fkbA RNA. This study uncovers a novel epigenetic RNAi-based epimutation mechanism controlling phenotypic plasticity, with possible implications for antimicrobial drug resistance and RNAi-regulatory mechanisms in fungi and other eukaryotes.

A single argonaute gene participates in exogenous and endogenous RNAi and controls cellular functions in the basal fungus Mucor circinelloides.

The mechanism of RNAi is well described in metazoans where it plays a role in diverse cellular functions. However, although different classes of endogenous small RNAs (esRNAs) have been identified in fungi, their biological roles are poorly described due, in part, to the lack of phenotype of mutants affected in the biogenesis of these esRNAs. Argonaute proteins are one of the key components of the RNAi pathways, in which different members of this protein family participate in the biogenesis of a wide repertoire of esRNAs molecules. Here we identified three argonaute genes of the fungus Mucor circinelloides and investigated their participation in exogenous and endogenous RNAi. We found that only one of the ago genes, ago-1, is involved in RNAi during vegetative growth and is required for both transgene-induced RNA silencing and the accumulation of distinct classes of esRNAs derived from exons (ex-siRNAs). Classes I and II ex-siRNAs bind to Ago-1 to control mRNA accumulation of the target protein coding genes. Class III ex-siRNAs do not specifically bind to Ago-1, but requires this protein for their production, revealing the complexity of the biogenesis pathways of ex-siRNAs. We also show that ago-1 is involved in the response to environmental signals, since vegetative development and autolysis induced by nutritional stress are affected in ago-1(-) M. circinelloides mutants. Our results demonstrate that a single Ago protein participates in the production of different classes of esRNAs that are generated through different pathways. They also highlight the role of ex-siRNAs in the regulation of endogenous genes in fungi and expand the range of biological functions modulated by RNAi.

A white collar 1-like protein mediates opposite regulatory functions in Mucor circinelloides.

Protein ubiquitylation plays a major role in the regulation of many cellular processes by altering the stability, localization or function of target proteins. CrgA is a protein of Mucor circinelloides that shows the characteristics of ubiquitin ligases and is involved in the regulation of carotenogenesis and asexual sporulation in this fungus. CrgA, which belongs to a poorly characterized group of proteins present in almost all eukaryotes, represses carotenogenesis through the proteolysis-independent mono- and di-ubiquitylation of Mcwc-1b, a White Collar-1-like protein which, when it is non-ubiquitylated, activates carotenogenesis. Using a proteomic approach, this work shows that the regulation of M. circinelloides vegetative development by CrgA is also mediated by Mcwc-1b, although, in this case, the non-ubiquitylated Mcwc-1b form acts as a repressor. High levels of a protein that contains a classical Rossmann-fold NAD(P)H/NAD(P)(+) binding domain for NAD(P)H binding and is similar to NmrA NADP(H) sensor-like proteins occur when Mcwc-1b is inactivated by ubiquitylation. A role for this protein in the regulation of sporulation is suggested because its over-expression suppresses the sporulation defect in a crgAΔ mutant. NmrA-like proteins are repressors that interact with GATA transcription factors and have been shown to be related to cell differentiation in Magnaporthe oryzae and Dictyostelium discoideum. This proteomic approach also revealed that CrgA regulates the carbon and energy metabolism and that Mcwc-1b is the main, but not the only, target of CrgA.

Malic enzyme activity is not the only bottleneck for lipid accumulation in the oleaginous fungus Mucor circinelloides.

Commercial interest in microbial lipids is increasing due to their potential use as feedstock for biodiesel production. The supply of NADPH generated by malic enzyme (ME; NADP+-dependent; EC 1.1.1.40) has been postulated as being the rate-limiting step for fatty acid biosynthesis in oleaginous fungi, based mainly on data from the zygomycete Mucor circinelloides studies. This fungus contains five genes that code for six different ME isoforms. One of these genes, malA, codes for the isoforms III and IV, which have previously been associated with lipid accumulation. Following a strategy of targeted integration of an engineered malA gene, a stable strain overexpressing malA and showing high ME activity has been obtained, demonstrating the feasibility of this strategy to overexpress genes of biotechnological interest in M. circinelloides. This is the first report showing the integration and overexpression of a gene in Zygomycetes. Unexpectedly, the genetically modified strain showed a lipid content similar to that of a prototrophic non-overexpressing control strain, suggesting that another limiting step in the fatty acid synthesis pathway may have been revealed as a consequence of the elimination of malic enzyme-based bottleneck. Otherwise, the fact that prototrophic strains showed at least a 2.5-fold increase in lipid accumulation in comparison with leucine auxotrophic strains suggests that a wild-type leucine biosynthetic pathway is required for lipid accumulation. Moreover, increasing concentrations of leucine in culture medium increased growth of auxotrophs but failed to increase lipid content, suggesting that the leucine synthesized by the fungus is the only leucine available for lipid biosynthesis. These results support previous data postulating leucine metabolism as one of the pathways involved in the generation of the acetyl-CoA required for fatty acid biosynthesis.

Linking genome annotation projects with genetic disorders using ontologies.

Genome sequencing projects generate vast amounts of data of a wide variety of types and complexities, and at a growing pace. Traditionally, the annotation of such sequences was difficult to share with other researchers. Despite the fact that this has improved with the development and application of biological ontologies, such annotation efforts remain isolated since the amount of information that can be used from other annotation projects is limited. In addition to this, they do not benefit from the translational information available for the genomic sequences. In this paper, we describe a system that supports genome annotation processes by providing useful information about orthologous genes and the genetic disorders which can be associated with a gene identified in a sequence. The seamless integration of such data will be facilitated by an ontological infrastructure which, following best practices in ontology engineering, will reuse existing biological ontologies like Sequence Ontology or Ontological Gene Orthology.

Molecular tools for carotenogenesis analysis in the zygomycete Mucor circinelloides.

The carotene producer fungus Mucor circinelloides is the zygomycete more amenable to genetic manipulations by using molecular tools. Since the initial development of an effective procedure of genetic transformation, more than two decades ago, the availability of new molecular approaches such as gene replacement techniques and gene expression inactivation by RNA silencing, in addition to the sequencing of its genome, has made Mucor a valuable organism for the study of a number of processes. Here we describe in detail the main techniques and methods currently used to manipulate M. circinelloides, including transformation, gene replacement, gene silencing, RNAi, and immunoprecipitation.

Two distinct RNA-dependent RNA polymerases are required for initiation and amplification of RNA silencing in the basal fungus Mucor circinelloides.

RNA-dependent RNA polymerases (RdRPs) play key roles in the RNA silencing pathway in a number of organisms. They have been involved in the production of double-stranded RNA (dsRNA) molecules that initiate the silencing mechanism as well as in the amplification of the silencing signal. The roles of RdRPs from fungi in these processes are poorly described compared with other eukaryotes. RNA silencing in the zygomycete Mucor circinelloides exhibits uncommon features, such as induction by self-replicative sense transgenes and an amplification process associated with two size classes of antisense small interfering RNAs (siRNAs). To investigate the function of fungal RdRP proteins in initiation and amplification of silencing we have cloned and characterized two different rdrp genes in M. circinelloides. Functional analysis of rdrp(-) disruption mutants indicates that rdrp-1 is essential for initiation of silencing by sense transgenes by producing antisense RNA transcripts derived from the transgene, but it is not necessary for amplification of the silencing signal, whereas rdrp-2 is required for efficient accumulation of the two different classes of secondary siRNAs regardless the nature of the trigger. Our results provide evidence for a functional diversification of M. circinelloides rdrp genes in different steps of the same RNA silencing pathway.

Sporangiospore size dimorphism is linked to virulence of Mucor circinelloides.

Mucor circinelloides is a zygomycete fungus and an emerging opportunistic pathogen in immunocompromised patients, especially transplant recipients and in some cases otherwise healthy individuals. We have discovered a novel example of size dimorphism linked to virulence. M. circinelloides is a heterothallic fungus: (+) sex allele encodes SexP and (-) sex allele SexM, both of which are HMG domain protein sex determinants. M. circinelloides f. lusitanicus (Mcl) (-) mating type isolates produce larger asexual sporangiospores that are more virulent in the wax moth host compared to (+) isolates that produce smaller less virulent sporangiospores. The larger sporangiospores germinate inside and lyse macrophages, whereas the smaller sporangiospores do not. sexMΔ mutants are sterile and still produce larger virulent sporangiospores, suggesting that either the sex locus is not involved in virulence/spore size or the sexP allele plays an inhibitory role. Phylogenetic analysis supports that at least three extant subspecies populate the M. circinelloides complex in nature: Mcl, M. circinelloides f. griseocyanus, and M. circinelloides f. circinelloides (Mcc). Mcc was found to be more prevalent among clinical Mucor isolates, and more virulent than Mcl in a diabetic murine model in contrast to the wax moth host. The M. circinelloides sex locus encodes an HMG domain protein (SexP for plus and SexM for minus mating types) flanked by genes encoding triose phosphate transporter (TPT) and RNA helicase homologs. The borders of the sex locus between the three subspecies differ: the Mcg sex locus includes the promoters of both the TPT and the RNA helicase genes, whereas the Mcl and Mcc sex locus includes only the TPT gene promoter. Mating between subspecies was restricted compared to mating within subspecies. These findings demonstrate that spore size dimorphism is linked to virulence of M. circinelloides species and that plasticity of the sex locus and adaptations in pathogenicity have occurred during speciation of the M. circinelloides complex.

Endogenous short RNAs generated by Dicer 2 and RNA-dependent RNA polymerase 1 regulate mRNAs in the basal fungus Mucor circinelloides.

Endogenous short RNAs (esRNAs) play diverse roles in eukaryotes and usually are produced from double-stranded RNA (dsRNA) by Dicer. esRNAs are grouped into different classes based on biogenesis and function but not all classes are present in all three eukaryotic kingdoms. The esRNA register of fungi is poorly described compared to other eukaryotes and it is not clear what esRNA classes are present in this kingdom and whether they regulate the expression of protein coding genes. However, evidence that some dicer mutant fungi display altered phenotypes suggests that esRNAs play an important role in fungi. Here, we show that the basal fungus Mucor circinelloides produces new classes of esRNAs that map to exons and regulate the expression of many protein coding genes. The largest class of these exonic-siRNAs (ex-siRNAs) are generated by RNA-dependent RNA Polymerase 1 (RdRP1) and dicer-like 2 (DCL2) and target the mRNAs of protein coding genes from which they were produced. Our results expand the range of esRNAs in eukaryotes and reveal a new role for esRNAs in fungi.